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Abstract. The corrections to the Curie temperature T c of a ferromagnetic 
film consisting of N layers are calculated for N 3> 1 for the model of D- 
component classical spin vectors in the limit D — > oo, which is exactly soluble 
and close to the spherical model. The present approach accounts, however, 
for the magnetic anisotropy playing the crucial role in the crossover from 3 
to 2 dimensions in magnetic films. In the spatially inhomogeneous case with 
free boundary conditions the D — oo model is nonequivalent to the standard 
spherical one and always leads to the diminishing otT c (N) relative to the bulk. 



The application of the spherical model EJ to spatially-inhomogeneous magnetic 
systems such as ferromagnetic films with free boundary conditions by Barber and 
Fisher Q has revealed the unphysical behaviour of the solution being the consequence 
of the global spin constraint. The dependence T C (N) for a ci-dimensional hypercubic 
lattice infinite in d! = d — 1 dimensions and having N layers in the dth dimension 
has been found to be for d > 4 a non- monotonous function with a maximum, i.e., 
T C (N) for N ^> 1 was larger than in the bulk. Other singular features of the spherical 
model were found by Abraham and Robert || by considering the problem of phase 
separation (i.e., the domain wall formation). 

Besides the numerous publications using the spherical model for inhomogeneous 
systems in its original form (see, e.g., [Q, |5|), there is a work by Costache, Mazilu und 
Mihalache Q in which the global spin constraint was replaced for a ferromagnetic 
film by separate constraints in each layer. Although this model is less convenient for 
analytical calculations, it was shown that for d > 4 the value of T C (N) monotonically 
increases to its bulk value T c (oo), as it should from the physical grounds. Earlier 
Knops 0| had proved that in a general inhomogeneous situation the spherical model 
with a spin constraint on each lattice site is equivalent to the D-component classical 
vector model by Stanley || in the limit D — > oo. The latter is not only more physically 
appealing than the original spherical model, but it also allows one to take into account 
the spin anisotropy |0j and to produce the l/D expansions |l^, O, |l2| . A convenient 
tool to handle the Z?-vector model is the classical spin diagram technique [|l2| [l3| . 

The possibility of considering anisotropic systems makes the analytically soluble 
D = oo model, which can also be called for simplicity the spherical model, rather 
attractive for applications. In Ref. |Q it was used to investigate the role of fluctuations 
in the phase transition from Bloch to linear domain walls in biaxial ferromagncts. 
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Anisotropy also plays a crucial role for ferromagnetic films in the actual case d = 3. 
For N ^ oo the system is infinite in only d! = 2 dimensions, and cannot sustain a long- 
range order in the case of a contunuous spin symmetry. Correspondingly, Barber and 
Fisher Q have found diverging corrections to T c for N 3> 1 for the standard spherical 
model. On the other hand, for a purely 2-dimensional system (N = 1) T c tends to 
zero only logarithmically slow with vanishing anisotropy. One can expect that in the 
quasi 3-dimensional case (N ^> 1) the characteristic anisotropy required to support 
the long-range order should be extremely small. The calculation of the T c -corrections 
for 3-dimensional ferromagnetic films, which strongly depend on the anisotropy, is the 
main purpose of this work. 

The hamiltonian of the anisotropic classical _D-vector model can be written in the 
form 

H = -- ^ Jij ( rn zi m z:j + 7] ^ m ai m a j J , (1) 

ij \ a=2 ) 

where is the normalized £>-component vector, |mj| = 1 and r\ < 1 is the dimen- 
sionless anisotropy factor. In the mean field approximation (MFA) the Curie tem- 
perature of this model is T c MFA = J /D, where Jo is the zero Fourier component of 
the exchange interaction. It is convenient to introduce the dimensionless temperature 
variable = T/T C MFA and the reduced correlation function (CF) of transverse (a > 2) 
spin components: Sy = D(m a im a j) , which are well-behaved in the limit D — ► oo. 
Using the diagram technique for classical spin systems |l2[ 0, one arrives in 
the limit D —> oo at the closed system of equations for the average magnetization 
m, = (m z i) and the CF Sij. These are the magnetization equation 

'"- (; < X! 'V/'" /- ( 2 ) 

3 

the Dyson equation for the correlation function 

sw — 0GiSai + i]Gi S XijSjv (3) 

3 

and the kinematic relation playing the role of the spin constraint on a lattice site i 



rn 
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1. (4) 



Here Ay = Jij/ Jo and Gi is the one-site spin average D(m a im a i) /9 renormalized by 
fluctuations, which should be eliminated from the equations. 

In the Ising case (77 = 0) the influence of the fluctuations of the transverse spin 
components disappear. Since, additionally, the longitudinal fluctuations dying out 
as 1/.D are not present in the equations above, the situation is in this case exactly 
described by the MFA. Equation (|^) has for 77 = the trivial solution su = OGi] then 
from constraint (Q) one gets Gi = (1 — mf)/6, and the elimination of Gi in (Q) leads 
to a closed equation for the magnetization. 

In the homogeneous case = m and Gi — G are constants, and the equation 
(j^) can be easily solved with the help of the Fourier transformation, which results in 

Here vq is the unit cell volume and Ak = Ju./ Jo- For the d-dimensional hypercubic 
lattices vo = aft, clq is the lattice spacing, 

1 d 

Ak = -y^cos(aofcj) ) (6) 

i=l 



2 



and the lattice integral P(X) has the following properties: 



l + X 2 /{2d), AT<1 

(l/7r)ln[8/(l-X)], 1-X<1, d = 2 

' ' ' W 3 -c 3 (l-Xy/*, 1-X«1, d = 3 

W 4 -c 4 (l-X)\n[c'J(l-X)}, d = 4, 



where VF 3 = 1.51639 and W 4 = 1.23947 are the Watson integrals, c 3 = (2/tt)(3/2) 3 / 2 
and c 4 = (2/ir) 2 . For d > 5 the leading terms of the expansion of P{X) about X = 1 
are non-singular. Since in the homogeneous case the sum in the right-hand side of (^) 
equals to to, it is satisfied only if to = (above 9 C ) or G — 1 (below C ). Then from 
equation (^) one gets the temperature-dependent magnetization: 

m = (i-e/e c ) 1 / 2 , e < e c = i/p( v ). (8) 

In the isotropic case (r/ = 1) for d > 3 the value of the transition temperature in 
the bulk 9 C reduces to the well-known result 9 C — 1/W [[[]. For d = 2 one obtains 
6 c {rj) = 7r/ ln[8/(l — 77)] vanishing for 77 — > 1. In the Ising case ?; = the MFA result 
C = 1 is reproduced. 

To solve the equations of the spherical model for a d-dimensional hypercubic 
ferromagnetic film it is convenient to use the Fourier representation in d! = d — 1 
translationally-invariant dimensions and the site representation in the dth dimension. 
The Dyson equation (|^) for the Fourier-transformed CF a nno (k) takes on the form of 
a system of the second order finite-difference equations 

2b n a n - er„ + i - <7„_i = [2d6/rj)8 nno , n = 1, 2, N, (9) 

where the mute index uq of a was dropped. For the free and periodic boundary 
conditions (fbc and pbc) in (JoJ) we set 

cto = o-jv+i = 0, (fbc) 

(T = <j n , <?n+i = 01, (pbc, N > 3). 

The coefficient b n in (^) reads 

b n = 1 + d^G,,)- 1 - 1] + (T(l - At), (11) 

where is given by @ with d ^ d' . The magnetization equation (||) takes on the 
form 

26„to„ - m„ + i - m„_i = (12) 

with b n = b n (i] = 1, k = 0) and the boundary conditions similar to (|l0|). The constraint 
equations (Q) can now be written as 

1 f d d 'k 

s nn + m 2 n = l, s nn = aq y ?2^p crnn(k), n = 1,2, ...,7V. (13) 

The solution of the equation (||) is governed by the effective k-dependent correlation 
length, which in the long wavelength region, a^k <C 1, is given by 



r c (k) = ao/^diivG)- 1 - 1] + (a fc) 2 , (14) 

and which should be compared with the film thickness L = Nag. In the region 
of parameters r c {k) <C L one can expect the d-dimensional quasi-bulk behaviour 
perturbed due to the finite L. In the opposite limit a behaviour corresponding to the 
reduced dimensionality d' = d — 1 is to be expected. For d > 3 in situations where 
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the finite-size corrections to 9 C are small, the main contribution to the integral (13]) 
comes from the region a^k ~ 1. For such wave vectors the correlation length (|l4) is 
of the order of the lattice spacing cto, and cr„„(k) are the functions of G n in several 
neighbouring layers. Then from the constraint equations (|l3|) it follows (at least in the 
paramagnetic state, m n — 0) that the inhomogeneity of G n in the fbc case is confined 
to the boundary regions n, N — n<n c ~ 1. Due to this inhomogeneity an analytical 
solution of the problem is possible only in limiting cases. In the subsequent we shall 
restrict ourselves to the calculation of the Curie temperature 9 C of ferromagnetic films 
with N > 1. 

In the Ising limit r] = we have G = 1/9 above 9 C , and 9 C can be found in the fbc 
case from the condition that the determinant of the linear system of equations ( |l2] ) 
turns to zero. This leads to the MFA result 

c = l-~(l-cos-^_y (15) 



d V N + 1, 

For the model with pbc G is also independent of n due to the symmetry of the 
problem, and for 9 < C , where m ^ 0, one finds G = 1 from (0) or (|l2|). The 
homogeneous solution of the finite-difference equation @ with b n = b has the form 
cr n = Ci/z" + C2fi~ n , and the result for the one-layer correlator a nn reads 



*nn(k)= f 1 | +/i - J V , M = 6+V / &^1. (16) 

r)Vb 2 -11- fJ, JV 

In the region 1 — r\ <C 1 and aofc <C 1 this expression has the limiting forms 
2d0 1 



-JV 



0"nn(k) 



r;iV 2d(l-7j) + (aofc) 2 ' 

L/r c (k) > 1 



L/r c (fc) = AV2d(l - '7) + (aofc) 2 « 1 



»7 v /2d(l-r/) + (a fc) 2 ' 

(17) 

demonstrating the crossover from d- to d'-dimensional behaviour mentioned above. 
The second of these limiting expressions corresponds to the bulk and can also be 
obtained by the integration of the bulk CF Sk (0) over the dth component of the wave 
vector. For N 1 and d — 3 the integral in ( |l3| ) with cr nn (k) fl6| ) can be calculated 
analytically, and the result for 9 C reads (pbc) 

0- 1 ^W 3 + ^-\n 1 ; (18) 

ttjV 1 - exp[-AV6(l - v)] 

In the limit of extremely small anisotropies l—rj the transition temperature 9 C becomes 
logarithmically small: 

but this limit is very difficult to reach for iV > 1. The minimal value of 1 — rj 
required to support 9 C ~ 1 diminishes exponentially fast with the increase of N: 
1 — rj* ~ jV -2 cxp(— 27tjV/3). For d = 4 the results have the form (pbc) 

^4+3^2, 1 « JV 2 « 1/(1 - rj) 

4\2(1 - n)} 1/4 , ^ 

W4+ (7rjv)3/ 2 exp[-2JVv/2(l-7?)], iV 2 (l-7?)»l. 

The first of these limiting expressions coinsides with that of Barber and Fisher j|. 
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Now we proceed to the investigation of the more complicated case of a ferromag- 
netic film with free boundary conditions. Here the solution of the Dyson equation 
with n = no can be represented by the recurrence formula 



2d9 1 1 



n 2b n — atn — a'n " + 2b n - a n ' 1 2b n - a' n 



(21) 



with the initial conditions 

ai = a' N = 0, a 2 = l/(2&i), a' N _ x = l/(2b N ). (22) 

Now all quantities G n entering b n ( |ll|) can be determined numerically as functions of 
9 from the N constraint equations (|13|). Finally, 6 C can be found from the condition 
Dn — 0, where Dn is the determinant of the linear system (|l2|). The problem can 
be solved analytically in two limiting cases depending on the value of N 2 (l — rf) (see 



(17)). In the limit 7V 2 (1 — 77) ^> 1 the system shows a <i-dimensional (bulk) behavior 
in the whole range of k, and in the main part of a sample all cr n „(k) are equal to 
each other and determined by the value of G far from the boundaries. Indeed, in 
this region the recurrence relations in (^) converge in the depth of the sample to 
a = a' = b — \Jb 2 — 1, which leads to the bulk expression for cr nra (k) analogous 
to the second one in @. The value of G in the depth of the film is in our limit 
insensitive to its behaviour in the boundary regions n, N — n<n c ~ 1 and can be 
found from the condition = using G n = G = const. To see that, one can simply 
cut the boundary regions and require _Djv_ Jlc = 0, which introduces corrections of 
the order 1/N <C 1. The calculation analogous to that in the MFA case yields G = 
l + (n/N) 2 /(2d). After integration over k in ( |l3| ) one arrives at the obvious expression 
Snn = 0GP(rjG) (see (H)), and the value of 9 C determined from the condition s nn = 1 
reads (fbc, N 2 (l _ n ) > 1) 



1 



1 



For d — 3 the limiting forms of I(r]) obtained from (Q) are given by 

l + n 2 /d, ?7<Cl 
I(r)) = < (3/2) 3 / 2 1 (24) 

— — , 1 - 77 < 1. 

irW 3 V^V 

One can see that (|23| ) generalizes the MFA result (^5|) , and for 1 — 77 <C 1 corrections 
to 9 C due to the finite-size effects arc much greater than in the MFA. 

For d > 5 the derivative P' (77) is finite for 77 — > 1 (see (^)) and the results obtained 
above can be applied for all values of n. The reason for this is that the region of small 
wave vectors, k<k^ = a^ 1 /N, where for A 2 (l — 77) < 1 the quasi-bulk expression for 
the CF er, m (k) becomes invalid (see (|l7|)), is suppressed by the phase- volume factor in 
the integral (|l3|). The marginal case is d = 4, where for 1 — 77 < (ao/c^v) 2 ~ 1/A^ 2 with 
the logarithmic accuracy it is sufficient to calculate the integral over the Brillouin 
zone down to A; at. As a result one gets (fbc) 



W* + ^hiN + o(J^y 1 « N 2 & 1/(1 - 77) 
Wi+ 2^ ln ^ N 2 (1- V )<1, 

2 



(25) 



(cf. (20)). An asymptotic dependence of the type hi(N)/N in the isotropic limit 



with a coefficient close to unity was obtained numerically for the spherical model with 
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the layer-constraint in ||. To the contrast, for the standard spherical model [Q with 
fbc 6- 1 = Wa + a/N with a < §. 

For a 3-dimensional ferromagnetic film in the limit 7V 2 (1 — rf) <C 1 the leading 
correction to the k-integral (|l^) and hence to 9 C comes from the longwavelength region 
k fc/v = o-o 1 /N, where a nn (k) behaves 2-dimensionally. The form of cr n „(k) in this 
region can be determined from the general formula ( pl| ) . Beyond the narrow boundary 
regions the quantities b n = b n (rj = l,k = 0), etc., satisfy b n = a n = a' n = 1, and 
the values of a n and a' n can be found from the recurrence relations (^TJ) with the 
help of the expansion with respect to small 1 — r/ and (aok) 2 . As a result on gets the 



same expression (17) in the same wave vector range k < fc^r . This region yields the 
contribution of the order (1/N) ln[l/(A 2 (l — ij))] into s nn (|l3|). The contribution of 
the region k>k^ into the correction to s nn can be estimated in the following way. 
For 1 - ij « 1 the finite-size correction described by (^) and (|24|) comes from the 
region of small wave vectors a^k ~ y/1 — r\. In our case, however, ciofc/v 3> \/l — ry, 
and the corresponding contribution is reduced to the value of the order 1/N. With 
the logarithmic accuracy the latter can be neglected in comparison to that of the 
2-dimensional region k<ltN- The final result for 9 C of the 3-dimensional model with 
free boundary conditions can be written as (fbc, A^ 2 (l — rj) <C 1) 

^^ w ^^mhv) +0 {^)- (26) 

The similarity of this result with ([l8]) is not surprising since in the relevant region, 
k'&kjy, where r c (k) > L, all N layers are strongly correlated with each other and the 
type of boundary conditions plays no role in the leading approximation. 

The author thanks Hartwig Schmidt for valuable discussions. The financial sup- 
port of Deutsche Forschungsgemeinschaft under contract Schm 398/5-1 is greatfully 
acknowledged. 
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